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I.   INTRODUCTION 

The  perishable  product  is  beginning  to  receive  attention 
in  the  inventory  problem.   Often,  it  has  a  fixed  lifetime 
after  which  it  becomes  useless  to  satisfy  demands.   News- 
papers, medical  drugs,  human  blood,  photographic  films  and 
some  foodstuffs  are  common  examples. 

An  inventory  problem  for  perishable  products  is  often 
treated  as  a  single-period  inventory  problem,  sometimes 
called  a  newsboy  problem.   The  newsboy  problem  was  introduced 
by  Morse  and  Kimball  in  1950  and  the  development  of  it  has 
been  presented  in  several  journals  and  publications.   A 
brief  history  of  it  was  given  by  Masuda  [Ref .  1] . 

A  well-known  characterization  of  the  newsboy  problem  is: 

The  newsboy  must  decide  how  many  papers  to  purchase 
for  resale  during  the  day.   If  he  buys  too  many  papers — 
that  is,  more  papers  than  people  will  purchase  from  him — 
he  will  have  papers  left  at  the  end  of  the  day  and  will 
incur  a  loss  on  each,  since  old  newspapers  have  little 
salvage  value.   If,  on  the  other  hand,  he  purchases  too 
few  papers  for  resale,  he  will  sell  out  early  and  his 
later  customers  will  not  be  able  to  buy  from  him;  he  thus 
incurs  a  cost  (lost  profit)  for  each  paper  he  could 
have  sold,  but  didn't.   The  newsboy's  problem  is  to  de- 
cide how  many  papers  to  have  on  hand  at  the  beginning  of 
the  day  so  as  to  minimize  the  losses  (due  to  surplus  or 
shortage)  at  the  end  of  the  day.   [Ref.  2] 

Finding  good  solutions  to  the  newsboy  problem  depends 
upon  the  amount  of  information  about  the  magnitude  of  demand 
X. 

In  recent  research  dealing  with  newsboy  ordering  deci- 
sions where  little  information  existed  about  demand,  it  was 


assumed  that  at  the  end  of  each  period  the  demand  for  that 
period  would  be  known  [Refs.  1,3,4].   In  order  for  this  to 
happen,  one  would  need  unlimited  supply  or  some  way  to 
measure  the  demand  from  unsatisfied  customers.   In  practice, 
this  is  often  not  the  case.   For  example,  when  the  newsstand 
is  empty,  the  later  customers  will  not  stop  at  that  news- 
stand.  Accordingly,  we  do  not  know  how  many  more  customers 
would  like  to  buy  that  newspaper.   Thus,  we  may  have  frequency 
data  on  demand  only  up  to  the  value  of  supply. 

This  thesis  is  intended  to  provide  methods  permitting 
solutions  to  the  newsboy  problem  when  we  only  have  demand 
frequencies  up  to  the  supply  value. 

In  Chapter  II,  two  procedures  will  be  proposed  to  esti- 
mate the  maximum  demand.   The  first  procedure  uses  the  Gompertz 
Curve  since  the  nature  of  the  demand  cumulative  distribution 
function  (c-d.f. )  can  be  viewed  as  a  growth  curve  with  the 
market  saturate   as  an  asymptote.   The  second  procedure  sug- 
gested is  the  four-parameter  Schmeiser-Deutsch  Distribution 
since  it  can  take  on  a  wide  variety  of  shapes  so  as  to  match 
many  distributions.   Also,  the  Schmeiser-Deutsch  Distribution 
was  found  to  be  a  suitable  procedure  for  the  case  of  limited 
supply. 

Using  the  estimation  of  maximum  demand  together  with  the 
unit  cost  of  surplus  and  the  unit  cost  of  outage,  the  estima- 
tion of  optimal  order  quantity  is  discussed  in  Chapter  III 
for  the  "under  uncertainty"  and  "under  risk"  cases. 
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Conclusions  and  applications  are  given  in  Chapter  IV. 
A  review  of  the  Schmeiser-Deutsch  Distribution  is  presented 
in  Appendix  A. 
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II.   ESTIMATION  OF  MAXIMUM  DEMAND 

Care  should  be  taken  to  estimate  the  maximum  demand  since 
the  optimal  ordering  quantity  for  uncertainty  depends  on  it. 
In  this  chapter,  two  alternatives  are  suggested  to  estimate 
the  maximum  demand  when  we  have  unlimited  supply.   Then,  we 
will  show  how  one  of  these  methods  may  be  used  to  estimate 
maximum  demand  when  there  is  no  data  on  unsatisfied  demand 
(the  case  of  limited  supply) . 

The  procedures  that  appear  suitable  and  realistic  are 
(1)  using  the  Gompertz  Curve  and  (2)  using  the  Schmeiser- 
Deutsch  Distribution.  We  will  discuss  the  Gompertz  Curve 
first. 

A.   USING  THE  GOMPERTZ  CURVE  TO  ESTIMATE  THE  MAXIMUM  DEMAND 

In  the  terminology  of  economics,  the  demand  c.d.f.  can 
be  considered  as  a  "growth  curve"  since  it  has  a  "market 
saturate"  as  an  upper  asymptote.   The  Gompertz  Curve  is  the 
most  generally  used  growth  curve.   Its  formula 


Y   =   kaD   ,  (1) 


where  X  is  the  time  series  and  Y  is  the  sales  amount,  de- 
pends on  the  values  of  the  parameters  k,  a  and  b  (or  their 
logarithms) ,  it  may  take  on  any  one  of  a  variety  of  shapes 
[Ref .  5]  . 
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+.   X 


Figure  1.   The  Gompertz  Curve  with  0<a<l,  0<b<l 

We  are  interested  in  using  the  Gompertz  Curve  with  0  <  a  <  1 
and  0  <  b  <  1  to  approximate  the  c.d.f.  of  demand.   This  form 
of  the  Gompertz  Curve  is  known  to  be  asymptotic  to  k.   When 
we  fit  the  Gompertz  Curve  to  the  observed  cumulative  fre- 
quency of  the  data,  k  should  theoretically  be  1.0.   In  prac- 
tice, however,  it  may  depart  slightly  from  1.0.   To  compensate 
for  this,  we  use  the  following  form  to  approximate  the  c.d.f.: 


Y      b 

k  =   a 


X 


The  following  description  is  taken  from  [Ref .  5] : 
The  parameters  k,  a  and  b  are  obtained  as  follows: 

1)  Be  certain  that  the  number  of  observations  is  divisible 
by  3,  i.e.,  there  are  3N  periods  of  base  data. 

2)  Take  the  common  logarithms  of  each  Y. 

3)  Add  the  first  N  log  Y's  to  obtain  £,  log  Y;  then  the 
second  N  data  points  to  obtain  E2  log  Y;  and  the  last  N 
points  to  obtain  l  _  log  Y. 

4)  Set  the  smallest  observation  to  be  unity. 
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5)  Substitute  in  the  following  formulas 


Z^log  Y  -  Z5log  Y 
bN   =    * ,,_  „ g^T=  »   ,  (2) 


Z  log  Y  -  Z  log  Y 


log  a  =   (Z9log  Y  -  Z .  log  Y)  ■  (b  -  1]     ,       (3) 
z         x        (b-l)z 


and 


log  k  =   i-  •  (Z  log  Y  -  (= — — )  •  log  a).     (4) 
N      L  b-1 


With  values  of  k,  a  and  b,  one  can  compute  Y  for  any 
value  of  X  using  the  formula  of  the  Gompertz  Curve.   Also, 
the  formula  can  be  derived  from  the  Gompertz  Curve  to 
generate  X  from  the  Gompertz  Curve  as: 


log(lo<?   I    -    lo<?   k) 

X   =  — 122* (5) 

log  b 


In  this  chapter,  we  wish  to  estimate  the  maximum  demand, 
If  the  demand  distribution  has  an  infinite  upper  tail,  we 
can  not  hope  to  estimate  the  100   percentile  of  demand. 

+■  h 

We  will  be  content,  therefore,  with  the  99    percentile  of 
demand.   This  is  found  by  setting  (Y/k)  =  .99  or  Y  =  .99k 
in  equation  (5),  i.e., 


.   ,-.0044, 
log  (  t ) 

Xm^  will  be   =   J-°%  *     ,  (6) 

max  log  b 
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which  may  be  considered  to  be  the  estimator  of  the  maximum 
demand  if  the  smallest  observed  X  is  unity.   (If  not,  the 
maximum  demand  should  be  adjusted  to  be  X  plus  the  smallest 
X  minus  one. ) 

As  an  example  assume  that  the  demands  for  newspapers 
have  been  10,  6,  9,  7,  5,  13,  11,  7,  8,  8.   Using  the  Gompertz 
Curve  to  estimate  the  maximum  demand,  we  can  preceded  as 
follows: 

Let  X  be  the  demand  for  newspapers,  and  Y  be  the  cumula- 
tive relative  frequency  of  X. 

TABLE  I 
DEMAND  AND  ITS  CUMULATIVE  RELATIVE  FREQUENCY 


5  0.1 

6  0.2 

7  0.4 

8  0.6 

9  0.7 

10  0.8 

11  0.9 

12  0.9 

13  1.0 


Since  Y  is  a  growth  curve,  it  can  be  approximated  by  the 

bX 
Gompertz  Curve,  Y  =  ka   .   To  determine  the  appropriate 

values  of  k,  a  and  b  using  the  method  described  above: 

1)   From  these  data,  there  are  9  distinct  values  between 

the  smallest  observation  (5)  and  the  largest  observation  (13), 

and  9  is  divisible  by  3.   (Here,  N  =  9/3  =  3.) 
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2)  In  the  third  column  of  Table  II,  Y's  are  converted  to 
log  Y's. 

3)  Values  of  Z ,  log  Y,  Z„log  Y  and  Z.log  Y  are  shown  in 
the  fourth  column  of  Table  II. 

4)  Since  the  smallest  observed  X  is  not  1.0,  the  maximum 
demand  should  be  adjusted  by  adding  the  difference  between  the 
smallest  observed  X  and  1.0.   This  is  because  the  Gompertz 
Curve  assumes  the  smallest  observation  is  unity. 

5)  From  equation  (2) ; 


bN   =   ((-. 0915)-(-. 4737) )/( (-. 4737) - (-2 . 0969) 


=   0.2355  . 
Therefore, 

b   =   0.6175  . 
From  equation  (3)  ; 

log  a   =   (-. 4737-(-2. 0969)  )  (.6175-l)/(. 2355-1) 
=   -1.0622  . 


Therefore, 


a  =   0.0867 


From  equation  (4)  ; 


log  k   =  i(-2. 0969-(. 2355-1) (-1.0622)/(. 675-1) ) 


=   0.0087. 
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Therefore, 

k   =   1.0203  . 

I.e.,  the  Gompertz  Curve  is  Y  =  (1. 0203) ( 0 . 0867) (0 ' 6175)  , 
Using  the  99th  percentile,  the  maximum  demand  is 


X      =   11.3964  +  (5-1)   =   15.3964 
max 

-   15. 


TABLE  II 
COMPUTATIONS  IN  FINDING  A  GOMPERTZ  CURVE 

ZLog  Y 


X 

Y 

Log    Y 

5 

0.1 

-1.0000 

6 

0.2 

-0.6990 

7 

0.4 

-0.3979 

8 

0.6 

-0.2218 

9 

0.7 

-0.1549 

10 

0.8 

-0.0969 

11 

0.9 

-0.0458 

12 

0.9 

-0.0458 

13 

1.0 

0 

E^log  Y  =  -2.0969 


Z2log  Y  =  -0.4737 


Z3log  Y  =  -0.0915 


A  program  for  the  TI-59  to  determine  the  maximum  demand 
when  the  smallest  observed  X  is  unity,  and  to  determine  the 
parameters  (k,  a  and  b) ,  is  given  in  Appendix  B  (Program  Bl) . 
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Experimentation  with  simulated  data  seems  to  indicate  that 
the  Gompertz  Curve  procedure  works  well  when  supply  is  un- 
limited, but  not  nearly  so  well  when  supply  is  limited. 

B.   USING  THE  SCHMEISER-DEUTSCH  DISTRIBUTION  TO  ESTIMATE  THE 
MAXIMUM  DEMAND 

K.  Pearson  pointed  out  that  a  theoretical  distribution  must 
have  at  least  four  free  parameters,  in  order  to  fit  the  four 
characteristics  (location,  dispersion,  skewness,  and  kurtosis) 
adequately  [Ref.  6].   Kottas  and  Lau  recommends  the  four- 
parameter  Schmeiser-Deutsch  Distribution  as  most  suitable 
[Ref.  6] .   The  details  of  the  Schmeiser-Deutsch  Distribution 
are  given  in  Appendix  A. 

The  Schmeiser-Deutsch  Distribution  has  a  finite  upper 
tail.   Furthermore,  if  X  is  a  random  variable  having  the 
Schmeiser-Deutsch  Distribution,  there  exists  a  closed-form 
expression  for  the  maximum  value  of  X.   This  expression  is 
a+b(l-d)  ,  where  a,  b,  c  and  d  are  the  four  free  parameters. 

Using  the  data  in  the  previous  example,  we  will  employ 
the  Schmeiser-Deutsch  Distribution  to  estimate  the  maximum 
demand.   We  first  estimate  the  parameters  using  the  algorithm 
in  Appendix  A: 

1)  The  mode  m  =  (7+8)/2  =7.5. 
Set  a  =  7.5. 

2)  Set  d  =  (Y(7)+Y(8) )/2  =  0.5. 

3)  Select  the  two  desired  quantiles  which  satisfy  equation 
(9) ,  namely, 
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p1   =   0.2  ,   x1   =   6 


and 


p2   =   0.9  ,   x2   =   11 

Substitute  in  equation  (10),  yielding  c  =  2.9453. 
4)   Substitute  in  equation  (11),  yielding  b  =  52.0120. 
Then,  the  maximum  demand  is 

X     =   a  +  b(l-d)c 
max 

=   14.2529   -   14. 

This  compares  very  well  with  the  result  obtained  by  using 

the  Gompertz  Curve  which  was  X    =15.39  64.   For  these  data, 
c  max 

there  is  no  substantial  difference  between  these  two  alterna- 
tive methods  of  estimating  X 

3   max 

A  program  for  the  TI-59  to  determine  the  maximum  demand 
and  the  four  parameters  of  the  Schmeiser-Deutsch  Distribution 
is  given  in  Appendix  B  (Program  B2). 

C.   ESTIMATING  MAXIMUM  DEMAND  WHEN  THERE  IS  NO  DATA  ON 
UNSATISFIED  DEMAND 

In  the  previous  work  it  was  implicitly  assumed  that  the 

value  of  demand  for  each  period  was  obtainable  even  if  demand 

is  greater  than  supply.   This  rarely  happens  in  the  real 

world  unless  backorders  are  taken.   Speaking  of  the  "newsboy 

problem",  the  newsboy  can  collect  the  data  in  the  range  of 
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demand  from  zero  to  ordering  quantity  (Q)  only.   When  the 
newspapers  are  sold  out,  he  will  not  know  how  many  more  were 
in  demand  for  that  day.   If  he  is  recording  the  sales,  his 
observed  data  on  that  day  will  be  X  =  Q  even  if  the  demand 
is  actually  X  >_  Q. 

The  Schmeiser-Deutsch  Distribution  can  be  easily  used  to 
investigate  the  case  when  data  for  demands  greater  than  the 
quantity  Q  are  not  available.   This  is  because  its  parameters 
can  be  estimated  using  only  the  mode  and  two  selected  distinct 
percentiles . 

As  an  example,  suppose  that  the  newsboy  orders  only  10 
newspapers  for  each  day.   His  sales  data  for  a  10  day  period 
are  10,  6,  9,  7,  5,  10,  10,  7,  8,  8.   (Instead  of  demand, 
which  was  10,  6,  9,  7,  5,  13,  11,  7,  8,  8.)   The  sales  data 
and  their  cumulative  relative  frequency  will  be: 

TABLE  III 
SALES  DATA  AND  THEIR  CUMULATIVE  RELATIVE  FREQUENCY 

X  Y(X) 

5  0.1 

6  0.2 

7  0.4 

8  0.6 

9  0.7 
10  1.0 
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Since  Pr(X  =  10)  is  actually  Pr  (X  >_  10),  we  will  ignore 
this  data  point. 

Following  the  algorithm  in  Appendix  A  for  the  rest  of  the 
data  points: 

1)  Set  a  =  7.5. 

2)  Set  d  =  0.5. 

3)  Select  p,  =  0.6,  x,  =  8,  p2  =  0.7,  x2  =  9  which 
satisfy  equation  (9) . 

Substitute  in  equation  (10),  yielding  c  =  1.5850. 

4)  Substitute  in  equation  (11),  yielding  b  =  19.2279. 

From  this,  the  maximum  demand  is  X    =13.9  09  3  which  is 

max 

close  to  X    when  the  newsboy  has  unlimited  supply, 
max 

With  this  X    ,  one  can  determine  Q*.   Procedures  to  do 
max 

this  are  discussed  in  the  next  chapter. 
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III.   ESTIMATION  OF  OPTIMAL  ORDER  QUANTITY 

In  the  previous  chapter  it  was  shown  how  to  estimate  the 

maximum  demand  (X    ) .   In  this  chapter,  the  estimate  of  the 

max  r 

optimal  order  quantity  (Q*)  for  the  case  of  uncertainty  and 
also  for  the  case  of  risk  will  be  discussed. 

A.   OPTIMAL  ORDER  QUANTITY  UNDER  UNCERTAINTY 

When  data  in  hand  is  not  sufficient  to  estimate  the 
possible  future  values  of  demand  or  when  the  distribution 
of  demand  is  not  known,  the  problem  of  deciding  on  Q  falls 
into  the  class  called  decisions  under  uncertainty. 

It  has  been  shown  that  the  optimal  decisions  under  uncer- 
tainty (Q*)  are  [Ref .  1] : 


i*   = 


c 

; )  •  X     ,   for  a  continuous 

c  +  c      max     ,     , 
s    o  demand  , 


and 


*        c 

Q    <   ( -r — )  (x    +1)   <   Q*+l  / 

c  +c    max 


for  a  discrete  demand,  where  c   is  the  unit  cost  of  surplus 
and  c   is  the  unit  cost  of  outage. 

Since  one  consequence  of  each  decision  is  another  obser- 
vation of  demand,  one  begins  under  uncertainty  and  uses  the 
uncertainty  optimal  order  quantity  until  demand  data  is  ade- 
quate to  estimate  the  distribution  of  demand  and  change  to 
the  optimum  order  quantity  under  risk. 
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B.   OPTIMAL  ORDER  QUANTITY  UNDER  RISK 

When  one  says  "under  risk",  one  means  that  the  probability 
distribution  of  demand  is  known  or,  at  least,  some  estimate 
of  the  probability  distribution  of  demand  is  available.   The 
latter  case  is  possible  using  the  Schmeiser-Deutsch  Distri- 
bution suggested  in  Appendix  A. 

The  optimal  order  quantity  under  risk  is  the  value  of  Q* 
such  that 


c 
F(Q*)   =   7 ,   for  a  continuous  demand 


c  +  c 


distribution , 


and 


c 
F(Q*-D   <      °c    <   F(Q*)  , 
s    o 

for  a  discrete  demand  distribution. 

It  should  be  noted  here  that  all  one  really  needs  to 
estimate  the  probability  distribution  of  demand  is  a  value 
for  the  [c  /(c  +c  )]    quantile.   The  inverse  c.d.f.  of  the 
Schmeiser-Deutsch  Distribution  can  easily  take  care  of  this 
problem. 

As  an  example,  suppose  that  the  newsboy  purchased  the 

newspapers  at  the  cost  of  IOC  each  and  he  can  sell  them  for 

2(K  each,  with  leftovers  having  a  salvage  value  of  5C  each. 

Let  p  be  the  ratio  [c  /(c  +c  )].   Therefore,  p  =  2/3. 

o   o   s  r 

Consider  the  data  from  the  previous  example  at  the  end 
of  Chapter  II.   These  data  are  adequate  for  a  decision  under 
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risk.   The  p   quantile  can  be  estimated  by  the  inverse 
c.d.f.  of  the  Schmeiser-Deutsch  distribution, 


F_1(p)   = 


a-b(d-p)     if   p  <  d 


a+b(p-d)     if   p  >  d  , 


where  the  four  parameters  were  found  to  be  a  =  7.5, 
b  =  19.2279,  c  =  1.5850  and  d  =  0.5. 

Therefore,  F   (p)  =  8.623,  which  is  the  value  of  Q*  for 
the  case  of  continuous  demand. 

For  the  case  of  discrete  demand,  the  optimal  order  quan- 
tity under  risk  is  the  value  of  Q*  such  that  F(Q*-1)  <  p  <  F(Q*), 
and  then  Q*  =  9  in  this  example. 

In  the  next  chapter  we  will  give  conclusions  and  applica- 
tions from  this  work. 
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IV.   CONCLUSIONS  AND  APPLICATIONS 

One  of  the  important  requirements  needed  to  obtain  the 
optimal  order  quantity  in  newsboy  problems  is  the  nature  of 
the  demand  distribution.   This  thesis  suggested  some  proce- 
dures which  can  be  used  to  estimate  the  demand  distribution 
even  if  data  on  unsatisfied  demand  are  not  available.   In 
the  case  where  it  is  hard  to  tell  what  common  distribution 
the  demand  should  belong  to,  the  Schmeiser-Deutsch  Distribu- 
tion was  felt  to  be  useful. 

It  is  fitting  at  this  point  to  explain  how  the  procedures 
from  this  and  other  research  on  the  newsboy  problems  would 
be  used. 

Suppose  we  are  starting  a  sequence  of  single-period 
inventory  ordering  decisions  but  have  no  information  on  the 
distribution  of  demand. 

One  way  of  employing  these  results  would  involve  three 
phases: 

1)   For  the  decisions  at  time  period  1  through  some  period 

m,  the  newsboy  would  probably  rely  upon  a  guess  at  the  value 

of  the  maximum  demand  (X    ) .   Considering  the  unit  cost  of 

max  * 

surplus  (c  )  and  the  unit  cost  of  outage  (c  ) ,  he  has  two 
alternatives : 

(a)  In  case  of  c   <<  c  ,  use  Q*  =  X 

s     o  max 

(b)  Otherwise,  use  the  uncertainty  decision  rule. 


Q*   =   [c  /(c  +c  )  ]  -X 

os    o    max 
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2)  For  time  periods  (m+1)  through  some  later  period  n, 

he  would  use  Q*  =  [c  /(c  +c  ) ] *X    with  the  value  of  X 

o   o    s    max  max 

obtained  from  the  Gompertz  Curve  or  the  Schmeiser-Deutsch 
Distribution,  using  data  from  the  first  m  periods.   Clearly, 
he  may  wish  to  revise  this  estimate  as  more  data  become 
available. 

3)  After  time  period  n,  he  can  switch  from  treating  the 

decision  as  a  decision  under  uncertainty  to  a  decision  under 

risk,  using  the  inverse  c.d.f.  of  the  Schmeiser-Deutsch  Dis- 

tribution  to  estimate  the  value  of  [c  /(c  +c  )]    quantile, 

o   o    s     ^ 

which  is  known  to  be  Q*  for  continuous  demand. 

It  should  be  noted  here  that  because  of  the  data  needed, 
m  must  provide  at  least  three  periods  where  he  did  not  sell 
out,  and  from  the  work  of  Yong-u  Sok,  n  should  provide  at 
least  five  periods  where  he  did  not  sell  out  [Ref .  4] . 

From  the  work  outlined  above  we  have  a  reasonable  means 
of  starting  a  sequence  of  newsboy  decisions  when  there  is  no 
information  about  the  demand  distribution  and  lost  sales 
are  unknown.   It  is  hoped  that  the  procedures  suggested  in 
this  paper  will  be  useful  to  those  interested  in  single- 
period  inventory  problems  and  their  solutions. 
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APPENDIX  A 
THE  SCHMEISER-DEUTSCH  DISTRIBUTION 

In  this  appendix,  the  Schmeiser-Deutsch  Distribution 
will  be  briefly  reviewed.   For  a  further  study  [Ref.  6]  and 
[Ref.  7]  are  suggested. 

Due  to  the  lack  of  convenient  methods  of  generating  random 
values  from  distributions  having  more  arbitrary  shapes ,  the 
Schmeiser-Deutsch  Distribution  was  developed  and  found  suit- 
able for  generating  random  values  from  many  common  statistical 
distributions . 

The  Schmeiser-Deutsch  Distribution  can  take  on  a  wide 
variety  of  shapes,  which  makes  it  useful  for  many  applica- 
tions.  Distributions  ranging  from  Bernoulli  trials  through 
U-shaped  distributions  to  the  uniform  distribution,  as  well 
as  heavier  tailed  distributions,  are  obtainable.   Also,  shapes 
ranging  from  symmetry  to  maximum  skewness,  as  measured  by  the 
third  moment,  are  possible.   The  Schmeiser-Deutsch  Distribu- 
tion can  have  values  for  the  first  four  moments  so  as  to 
match  any  of  the  well-known  distributions,  and  straight- 
forward parameter  determination  techniques  exist. 

A  most  favorable  characteristic  of  the  Schmeiser-Deutsch 
Distributions  is  the  existence  of  simple  closed-form  expres- 
sions for  both  the  c.d.f.  and  the  inverse  c.d.f.,  being 
respectively, 
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F(x)   =   p   = 


and 


d-[ (a-x)/b]1/c    if   (a-bdc)  <  x  <  a 


d+[  (x-a)/b]1//c    if   a  <  x  <  a+b(l-d)c 


(7) 


x  =   F_1(p) 


a-b(d-p) 


a+b(p-d) 


if    p  <  d 


if   p  >  d 


(8 


where  a,  b,  c  and  d  are  the  four  parameters  of  the  distribu- 
tion, and  p  is  a  percentile  [Ref.  7]. 

"The  location  and  the  spread  of  the  distribution  are 
determined  by  a  and  b,  respectively.   The  shape  of  the  dis- 
tribution, often  measured  in  terms  of  skewness  and  tail- 
weight  (or  peakness) ,  is  determined  by  -(d-p)   in  case  of 
p   d,  or  (p-d)   in  case  of  p  >  d.   This  factor  is  a  trans- 
formation from  the  rectangular  shape  of  the  uniform  distribu- 
tion to  the  shape  of  interest.   The  d  =  0.5  yields  a  symmetric 
distribution.   If  c  >  1,  then  skew  is  to  the  right  for 
d  <  0.5  and  to  the  left  for  d  >  0.5.  (If  c  <  1,  the  direction 
of  skew  is  reversed.)"  [Ref.  7] 

"The  distribution  possesses  two  characteristics  which  may 
adversely  affect  its  usefulness  for  specific  modeling  pur- 
poses.  The  first  characteristic  is  that  the  value  of  the 
density  function  evaluated  at  the  mode,  f(x  =  a),  assumes 
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only  three  values:   zero  for  0    c  <  1,  one  for  c  =  1, 
and  infinity  for  c  >  1.   The  second  characteristic  is  the 
truncated  tails  of  the  distribution."  [Ref.  7] 

"Despite  the  truncated  tails,  the  distribution's  versa- 
tility in  assuming  a  wide  variety  of  shapes  makes  it  a 
reasonable  model  for  a  wide  range  of  processes.   This  is  due 
to  the  ability  to  assume  the  essential  features  of  a  process, 
which  often  makes  the  truncation  insignificant. "  [Ref.  7] 

Parameter  Determination  [Ref.  7] 

The  mode  and  quantiles  are  used  in  estimating  parameters 
for  the  Schmeiser-Deutsch  Distribution.   The  Schmeiser-Deutsch 
Distribution  is  useful  in  two  common  situations.   The  first 
is  when  full  information  is  not  available,  but  such  quantities 
as  the  "most  likely  value"  and  the  "minimum"  and  "maximum" 
values  can  be  estimated.   The  second  situation  is  when  data 
is  in  histogram  form  where  F(m)  is  the  percentile  of  the 
mode  m.   The  two  desired  quantiles,  p,  =  F(x, )  and  p_  =  F(x2) , 
must  satisfy  either  of  the  following  two  expressions: 


( | a-x, |  <  |x2-a|   and   |d-p, |  <  |p2-d|) 

(9) 
( | a-x, |  <  | x2~a |   and   | d-p1 |  <  | p2~d | ) 


Although  computations  to  determine  the  four  parameters 
can  be  easily  handled  by  computer,  much  attention  has  been 
given  in  this  thesis  to  develop  "manual  procedures",  meaning 
manual  computational  capability.   The  TI-59  handheld  calculator 
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(characterized  by  low  price  and  surprising  capability  and 
compactness)  can  be  considered  as  a  "manual"  tool,  and  a 
program  for  the  TI-59  to  determine  the  four  parameters  is 
given  in  Appendix  B. 

An  algorithm  for  estimating  the  parameters  is  as  follows: 
[Ref.  7] 

1)  Set  a  =  m 

2)  Set  d  =  F(m) 

3)  Set 


=   In 


a-x. 


a-x. 


/  In 


d-p. 


d-p. 


(10) 


4)   Set  b  =  | a-x  |/|d-p. 


(11) 


5)   Due  to  the  existence  of  the  inverse  c.d.f.  in  closed- 
form,  the  value  of  X  corresponding  to  a  percentile  p  of  the 
distribution  may  be  found  by  inserting  p  into  equation  (8). 

Furthermore,  the  range  of  X  can  be  calculated  easily 

by  substituting  p  =  0  and  p  =  1  at  Step  5) ,  or  obtained 

directly  from  equation  (7)  since  X    =  a-bd   and  X    =  a+b(l-d) 

mm  max 

At  this  point,  one  should  determine  whether  or  not  the  ob- 
tained distribution  is  an  adequate  model  by  plotting  the 
p.d.f.  or  c.d.f.  [Ref.  7]. 
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APPENDIX  B 
TI-59  Programs 

This  appendix  presents  two  TI-59  programs  to  determine 
the  maximum  demand  and  the  parameters.   The  first  program 
is  using  the  Gompertz  Curve  with  the  restriction  that  the 
smallest  X  is  unity.   The  second  program  is  using  the 
Schmeiser-Deutsch  Distribution  with  the  two  desired  quantiles 
satisfying  equation  (9). 

Program  B, 

1.  Program  Description 

The  program  will  estimate  the  maximum  demand  using 

the  Gompertz  Curve,  and  will  also  determine  the  parameters 

bX 

(k,  a  and  b)  in  Y  =  ka 

2 .  User  Instruction 

Step    Procedure  Enter    Press      Display 

1  Enter  n  n        GO         n 

2  Enter  Y1,...,Y3n       Y1      ^751       Y± 


Y 
n 

IR/SI 

Z,log  Y 

Yn+1 

IR/SI 

Yn+1 

Y2n     ITT751        Z2log  Y 
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Step    Procedure  Enter      Press      Display 

Y2n+1        QSD         Y2n+1 


3n 


^751       23log  Y 


3  Compute  X  DEI         X 

v  max  ^         max 

4  Review  parameters 
(optional) 


a 

b 


The  maximum  demand  shown  on  the  display  is  for  the  case 
when  the  smallest  observed  X  is  unity.   Otherwise,  the  maxi- 
mum demand  should  be  adjusted  by  adding  X    -1,  where  X 

J  J  3   mm  mm 

is  the  smallest  observed  X. 


E3 

k 

0Q 

a 

ran 

b 
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Program  B2 

1 .   Program  Description 

Estimate  the  maximum  demand  using  the  Schmeiser- 
Deutsch  Distribution.   Also  determine  the  parameters  (a,  b, 
c  and  d)  in 


F(x)   = 


d  -  [— r— ]  /      if    a-bd    x  <_  a 


d  +  [^]1/C     if    a  <  x  <  a+b(l-d)c 
b  —   — 


2.      User   Instruction 


Step 

Procedure 

Enter 

1 

Enter  mode  m 

m 

2 

Enter  F (m) 

F(m) 

3 

Enter  p. 

Pl 

4 

Enter  x. 

Xl 

5 

Enter  p2 

p2 

6 

Enter  x? 

x2 

7 

Compute  X 
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